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Abstract 

The LSW theory of Ostwald ripening concerns the time evolution of the size 
distribution of a dilute system of particles that evolve by diffusional mass transfer 
with a common mean field. We prove global existence, uniqueness and continuous 
dependence on initial data for measure-valued solutions with compact support in 
particle size. These results are established with respect to a natural topology on 
the space of size distributions, one given by the Wasserstein metric which measures 
the smallest maximum volume change required to rearrange one distribution into 
another. 


1 Introduction 


The classical theory of Ostwald ripening, formulated by Lifshitz and Slyozov ||] and 
Wagner [5| concerns the evolution of the size distribution of a large number of small 
particles of one phase embedded in a matrix of another phase. Particles are assumed 
to be widely separated spheres that evolve by diffusional mass transfer with a common 
mean field. In the late stages of the phase transformation, diffusion is quasi-steady and 
the particle growth rate is determined by the mass flux at the particle boundary. The 
mass flux is proportional to the gradient of a potential that is harmonic, is proportional 
to curvature on the particle boundaries, and is close to constant in the mean field between 
particles. 

In appropriate units, it is found that any particle radius R(t) evolves according to 


dR 

dt 


V(R,Rc(t)): 


a 

Jr 



(i) 


where a is a constant and the critical radius R c {t ) is the same for all particles. The value 
of R c (t) is determined from conservation of mass. If mass changes in the diffusion field 
can be neglected, the particle volume is conserved and one finds that the critical radius 
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equals the average radius of currently existing particles. Particles with radius larger than 
R c (t ) are growing, and particles with smaller radius shrink and can disappear in finite 
time. 

Classically, the size distribution of particles is described by a particle radius distribu¬ 
tion n(t, R). This is a normalized number density that we may scale so that f^n(t, r ) dr 
is the number of (currently existing) particles with radius less than R, divided by the 
number N of initially existing particles. The number of particles with size between 
R.\(t) and Ri{t) for any two solutions of © is conserved, so n(t, R) should satisfy the 
conservation law 

d t n + d R (Vn) = 0, (2) 

where the critical radius is given by 

p oo / poo 

R c {t) = J Rn(t,R)dR / J n(t,R)dR. (3) 

The initial number density na(R) = n(0, R) satisfies Jq°° no(R) dR = 1 in this normaliza¬ 
tion. 

Our aim in this paper is to develop a satisfactory theory of well-posedness for the 
initial value problem for the particle size distribution. From the physical point of view, it 
is reasonable to suppose that a positive fraction of the particles can have the same radius, 
in which case the size distribution contains one or more Dirac deltas. Mathematically, 
the ideal is to allow the initial data no(R) dR to be an arbitrary probability measure such 
that the total volume f ( ^ ~7rR 3 no(R) dR is finite. 

It will be convenient to work with particle volume v instead of radius R, and to work 
with a cumulative number distribution function ip instead of the number density n. We 
say that 

p is the fraction of (initially existing) particles with volume > v. (4) 

As a function of volume v at time t. ip(t, v ) is a monotonically decreasing function which 
is left continuous at jumps with ip(t, 0) = 1, and J 0 °° tp(t, v) dv (the total volume) is 
independent of time. The particle volume distribution, defined by f(t, v ) dv = —dip(t, v ) 
for each fixed t, is formally related to n via f(t, v ) dv = n(t, R) dR. 

We normalize the time scale by the factor 47ra and let 9{t ) = (47ri2 c (t) 3 /3) -1 / 3 , so 
that the volume v{t) of any existing particle should satisfy 

^=A(M(f)):=« 1/3 0(i)-l. (5) 

If v(t) is a positive solution of (j|) on some time interval, then ip(t,v(t)) should remain 
constant. This means (p(t, v) should be a solution of the hyperbolic equation 

d t tp + A(v,6(t))d v <p = 0, (6) 

whose characteristics satisfy (|5|). The value of 0(t) is obtained from ip in terms of 
Riemann-Stieltjes integrals by 

poo j poo 

8(t) = J dip(t,v) / J u 1 / 3 dip(t, v) . (7) 
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The numerator is —1 times the quantity po(t) := lim^o v)> which is the fraction of 
initially existing particles that still exist at time t. 

It turns out to be still better to regard the volume v as a function of the fraction p, 
0 < p < 1. We take the map p 1 —► v(t,p) to be right continuous and decreasing with 
v(t , 1) = 0. Mathematically, given p(t,v) we obtain v(t,p) via the prescription 

v(t,x) = sup{y | p(t,y) > 2 ;} for 0 < x < 1 = max ip. (8) 

This is most easily understood when the size distribution corresponds to a finite number 
of particles. If we list the particle volumes in decreasing order, vq( t) > ... > vn-i (t), 
then v(t,ip) = Vj for p E \j/N, (j + 1)/N). We shall call <p 1 —> v(t,p) a volume ordering 
for the system at time t. 

For technical simplicity we shall assume that the particle volumes in the system 
are bounded. This seems reasonable physically, and corresponds to assuming that the 
particle volume distribution has compact support in v. We then introduce function spaces 
as follows. Let rcd{[ 0,1]) be the set of functions v: [0,1] —> M that are right continuous, 
decreasing, and satisfy u(l) = 0. (To be precise, we say v is decreasing if v(x\) < v(x 2 ) 
whenever x\ > X 2 , and similarly for increasing. A decreasing function need not be 
strictly decreasing.) The set rcd([ 0,1]) is contained in the space bdd([ 0,1]) of real-valued 
bounded functions on [0,1], equipped with the sup norm ||u|| = sup^, |u(</?)|. rcd([ 0,1]) is 
a complete metric space in the induced topology. 

If X is a Banach space and I C M is an interval, then C(I, X ) is the space of continuous 
X-valued functions on /, and L^ C (I) is the space of equivalence classes of measurable 
functions locally bounded on I, where two functions are considered equivalent if they 
agree almost everywhere. 

Our main results are the following. 

Theorem 1.1 (Global existence and uniqueness) Let vq E rcd([0.1]). Then there exist 
unique functions 9 E L^ c (0, 00 ) and v E C([0, 00 ), rcd([0, 1])), such that 

/ v(t,<p)d<p= / v 0 (<p)d<p 
Jo Jo 

for all t > 0, and 

v(t, p) = v 0 (p) + f {v(s, p) 1/3 9(s) - 1) ds 
Jo 

for all (t,p) such that v(t,p) > 0. 

Theorem 1.2 (Continuous dependence on initial data) Given positive constants T and 
Cq, there exists a positive constant C such that, if (vi,6\) and (^ 2 ,^ 2 ) are two solutions 
with the properties stated for ( v , 9) in Theorem 0 and */max(ui(0,0), ^ 2 ( 0 ,0)) < Co, 
then 

sup \\vi{t, •) - v 2 (t, -)|| < C'||ui(0, 0 - ^ 2 ( 0 , Oil- 

0 <t<T 

Consequently, the map vq v is locally Lipschitz from rcd([0,1]) into C([0, T],rcd([ 0,1])). 
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Our strategy to prove these results at the same time justifies a method of nu¬ 
merical approximation for the problem that has a direct physical interpretation. We 
first consider solutions that are piecewise constant, taking a finite number of values 
vo(t) > ... > VN-i(t), as is the case for a finite number of particles. We show that 
these solutions are determined on a succession of time intervals by solving finite systems 
of coupled ordinary differential equations with a number of components that decreases 


as the smallest particles vanish. Once we prove the continuity estimate in Theorem 1.2 


(at first for initial data near to each other), uniqueness is immediate and existence for 
general initial data in rcd([ 0,1]) follows by an approximation argument. 

The solutions constructed in Theorem |1.1| correspond to measure-valued weak solu¬ 
tions of the evolution equation 


d t f + d v (A(v,9(t))f) = 0 


(9) 


for the particle volume distribution. This means that at each time t, the formal expression 
f(t, v ) dv corresponds to a probability measure v t having compact support in [0, oo), the 
set of volumes. The notion of distance used in Theorem |1.2| has an interpretation as a 
natural metric on the space Vo of such probability measures. This metric measures the 
smallest ‘maximum volume change’ required to rearrange one volume distribution into 
another. Mathematically it is the L°° Wasserstein metric [|], [Ij, which we denote by d^. 
In section 3 we shall establish the relationship between v(t,cp) and z/j, and deduce the 
following result as a corollary of Theorems 0 and |1.2|. 


Theorem 1.3 Let Vo denote the set of probability measures on [0, oo) of compact support, 
with topology given by d^, the L°° Wasserstein metric. Given uq € Vo, there exists a 
unique 9 £ L“ c (0, oo) and a unique map t > v t that is locally Lipschitz from [0, oo) into 
Vo, such that (9, v) is a volume-conserving weak solution of ||j, in the sense that 

/•oo /•oo 

/ vdv t (v) = / vdv 0 (v) 

Jo Jo 

/•oo /*oo 

/ / d t ((t,v) + A(v,9(t))d v ((t,v) dv t (v)dt = 0 

io Jo 


for all t > 0, and 


for all smooth (: (0, oo) x (0, oo) —* R with compact support. 

Furthermore, given any T > 0, Co > 0, there exists C > 0 such that, if two such 
weak solutions (9i,o 1 ' 1 ^), (92, are given, such that the supports of and are 
contained in [0, Co], then 

sup doo^,^) < C'd 00 (^ 1) ,^ 2) ). 

0 <t<T 


It is arguably natural from the physical point of view to measure distance between 
volume distributions by using the Wasserstein distance as is done here. A physically 
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reasonable notion of distance should reflect in a plausible way the effect of small pertur¬ 
bations of the system on size distributions. In late-stage Ostwald ripening one imagines 
that the nucleation or destruction of large particles is unlikely. Thus the topology should 
not make it ‘easy’ to change the number of large particles. It is plausible, rather, that 
small perturbations to the system would involve small changes to particle volumes. These 
notions are captured here by the use of the sup norm distance between volume orderings, 
and this is equivalent to using the L°° Wasserstein metric to compare volume distribu¬ 
tions. 

In section 4 we briefly treat a related, but simpler, case that arises in LSW theory, in 
which mass variations in the diffusion field are not neglected. In this case it is not total 
particle volume that is conserved in time, but rather a quantity of the form 


iO(t) + f v(t,ip)d(p, 

J o 


where a > 0 is constant. The evolution of particle volumes is still given by (I) , but 9 is 
now determined directly from the conserved quantity. 

2 A priori estimates and well-posedness 


In order to prove the a priori estimate stated in Theorem |1.2| , we need a pair of lemmas 
that yield strengthened variants of Gronwall’s inequality. 

Lemma 2.1 Suppose G: [0, T] —> M is increasing with G(0) = 0, K > 0 is a constant 
and f: [0, T] —> M is continuous and satisfies 


0 < f{t) < K + / f{s) dG{s ), 0 < t < T. 

J o+ 

Then f(t ) < Ke G ® forO<t<T. 

Proof: Let 

U(t) = K+ [ f(s) dG(s), 

J o+ 

then U( 0) = I\ and U is increasing. To prove the lemma it suffices to show that e~ G U < 
K. Let {tj}’j = o be a partition of [0,T] and define 

A t= sup (tj - tj-fi), e(A t) = sup |/(f)-/0)|. 
l<j<n \t— s|<At 

Put Uj = U(tj ), Gj = G(tj). Then 

g-Gj+i Uj+i _ e -Gj jj. = e - G i+i(u j+1 - Uj ) - Uj(e~ G i - e ~ Gj+1 ) 


= e 


Gj+1 (/ J+1 ~ Uj { e ° j+1 ~ Gj ~ 


< e~ Gj+1 (f(tj) + e(A t) - Uj) (G j+ i - Gj) 

< e(A t)(G j+1 -Gj), 
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where we used that e x — 1 > x for all x and f(tj) < Uj. Summing, we find that 
e~ G iUj < K + e(At)Gj for all j. Since the partition is arbitrary, e(A t) can be made 
arbitrarily small and the result follows. 


Lemma 2.2 Suppose G: [0, T] —» R is increasing, and f: [0, T] —> M is continuous and 
nonnegative and increasing. Then as long as 0 < t + f(t) < T we have 


rt r f( 0 ) 

(G(s + f(s)) — G(s)) ds < 


10 


(G(f(0))-G(s))ds+ [* f(s)dG(s) 
J o 


where G(s) = G(s + f(s)). 


Proof: Let Q denote the quantity on the right hand side of the desired inequality. 
Observe that since G is increasing, we have that 


Q + / (G(s + f(s)) - G(s))ds >Q+ G(s) ds - G(t + f(t))f(t). 

Jt Jt 

Since G(t + f(t)) = G(t), after integrating by parts in the Riemann-Stieltjes integral and 
cancelling boundary terms we find that the last right hand side equals 


*+/(*) 


G _ /' t +/W _ 

/ G{s)ds- / G{s)df(s)+ / G(s) 

Jo Jo Jt 

rf(o) rt _ ft+nv _ 

= — G(s) ds — / G(s) d(s + f (s)) + / ^( s ) 

Jo Jo Jo 

ft+f{t) 

= / (G(s + f(s)) — G(s)) ds. 

■Jo 


ft+f(t) 


ds 


ds 


Cancelling the part of the integral from t to t + f(t) finishes the proof. 

Next we establish some basic properties of solutions of the initial value problem as 
described in Theorem |l.2| . Fixing T > 0, we shall consider t € [0, T]. Let 9 € L°°(0,T) 
be positive and let v £ C([0,T],rcd([0, 1])) be such that 


for all t and 


/ v(t,<p)d<p= / v(0,<p)d<p 

Jo Jo 


i(t, tp) = v(0, ip) + / (u(s, (/?) 1 / 3 0(s) - 1) ds 


( 10 ) 

(ID 


whenever v(t , tp) > 0. By scaling, we may assume v(t, tp) dip = 1 for all t £ [0, T], 
From ( P ) it follows that t > v(t, ip) is Lipschitz and satisfies 


( 12 ) 
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for almost every t in any interval where v > 0. Since v 1 t 3 Q — 1 < — ^ for v < £o where 

e^" 1 = 8 ess sup 0(t), 
o <t<T 

it follows easily that if v{to,p) = 0 then v(t,p) = 0 for all t > to- 

We define v(t) = v(t, 0) = max,, v(t, p) and with the notation a A b = min(a, b) we 
define 

t{p) = inf{t E [0, T] | v(t, p) = 0} A T, 

(pit) = sup{y? E [0,1] | v(t,p) > 0}. 

The functions t and p are decreasing functions, and (p(t) > 0 for all t, since v(t, •) can 
never vanish identically by volume conservation. We call t(p) the vanishing time for 
v(t, p) at p if t{p) < T (but note that t(p) = T if v(T, p) > 0). 


Lemma 2.3 For almost every t E [0, T] we have 

0 < Q(t) = 1 ^ -< v{t) 2/3 < (e t u(0)) 2 ^ 3 . 

[ v(t,p) 1/3 dp 
J o 

Proof: Evaluate ( [Tl] ) at min(t, t(p)) and integrate over p E [0,1]. Changing the order 
of integration and using the fact that v(t(p), p) = 0 if t{p) < t, we obtain 


0 = 


rCp(t) r 1 

/ v(t,p)dp— / v(0,p)dp 
Jo Jo 


<P(t) 

1 /»min (t,i(cp)) 


(v(s, p) 1 J 3 0(s) — 1) ds dp 

ft r<p(s) 

/ / (v(s, p) 1 / 3 0(s) — 1) dp ds. 

Jo Jo 


io Jo 

ft rip(s) 


Since t is arbitrary the formula for d(t) follows. To get the inequalities, we use that 
p{t) < 1 and v 1 / 3 dp > h(t)~ 2 / 3 vdp. Then since dv/dt < v l i 3 9 < v we find 

v(t) < e t v{0). (13) 


Lemma 2.4 Whenever v(t\,p) < £q, we have dv/dt < — | for almost every t E [t\,t(p)\ 
and 

\ii(P) ~t)< v{t, p) <£o~ g (t ~ h) 

for all t E [ti, t ((/?)]. 
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Proof: v < £q implies u 1 / 3 # — 1 < — ^ almost everywhere, and the results follow easily. 


Corollary 2.5 There is a constant C = C(T, Cq) such that 

rt(v>) 


pm 

/ v(t, dt < C 

Jo 


for all p £ [0,1]. Furthermore, the function f3 given by 

!•&(*) 


mo 

(3{t) = / v(t,p)~ 2/3 dip 

Jo 


is finite for a.e. t £ [0, T) and J 0 T f3(t) dt < C. 


Proof: The first assertion follows directly from the estimates of the preceding lemma. 
The second follows from Fubini’s theorem. 


Our plan now is to first prove a restricted version of Theorem |l.2| , for two solutions 
that are initially close together. This restricted result will suffice to establish the existence 
and uniqueness theorem, after which the results of Theorem L2 without restriction can 
be proved. 

Proposition 2.6 Given T > 0, Cq > 0, there exist C > 0 and 5 > 0 such that the bound 
asserted in Theorem Li holds under the additional assumption that 

IMO, ■) — ^ 2 (0, -)II < 5 . 


To start the proof of this restricted version of Theorem 1.2, we suppose that T, Cq > 0 
are given and put 

£1 = (8e T Co) _1 . 

We suppose that (9i,vi) and (# 2 , V 2 ) £ L°°(0, T) x C([0,T\,rcd([0, 1])) are two solutions 
of © and ( 0 ) such that max(ui(0, 0), ^ 2 ( 0 , 0)) < Cq. We define 

M(t)= sup ||ui(s,-) -u 2 (s,-)|| 

0 <s<t 

and assume that M( 0) < e\. 

Lemma 2.7 There is a constant C\ = Ci(T, Co) such that for 0 < t < T we have 

M(t) < Ci (m(0) + J |0i(s) - 0 2 («)| ds ' 
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Proof: Fix p G [0,1]. We suppose that ti(p) > t 2 (p) without loss of generality. For 
t G [0,i 2 (v)} we nmy write 

vi(t,p)-v 2 (t,p) = vi(0,tp)-v 2 (0,ip)+ [ v 2 (s,p) 1/3 (0i(s) - 0 2 (s))ds 

Jo 

+ [ 9i(s)(vi(s,<p) 1/3 - v 2 (s,p) 1/3 )ds. 

Jo 

Using the bounds above for 9\ and v 2 , and the fact that |a — 6| < |a 3 — & 3 |/a 2 whenever 
a, b > 0, with U* = (e T Co) 1 ^ 3 we obtain the estimate 

\vi{t, ip) - v 2 (t, p)\ < \vi(0, ip)-v 2 (0,ip)\+C* [ \6i(s) - 9 2 (s)\ds 

Jo 

+Cl ( vi(s,^) _2/3 |?;i(s,^) -v 2 {s,ip)\ds. (14) 

Jo 

For t G [t 2 (ip),ti (</?)], we have v 2 (t, p) = 0 and may write 

vi(t, p) < v 1 (t 2 (p),p) + Cl [ vi(s,p) 1/3 dp. 

Jt2(ip) 

Using ( pD with t = t 2 (p ) to estimate v\(t 2 (p),p), we find that (P) is valid for all 
t G [0, ti (</?)]. Gronwall’s inequality then yields that 


exp 


—C* j vi(s,p) 2/3 ds' S j \v!(t,p) -v 2 (t,p)\ 

< M0,y>) -v 2 (Q,p)\ + C* [ |0i(s) - 6» 2 (s)| ds. 

Jo 

Using Corollary |2.5| completes the proof. 


Lemma 2.8 Suppose M(t) < E\ for 0 < t < r. Then 

\<Pi(t) - p 2 (t )| < pi(t) - pi(t + 2 M(t)) + p 2 (t) - <p 2 (t + 2 M(t)). 
as long as t + 2 M(t) < r. 

Proof: Fixing t. by relabeling we can assume <pi(t) < <p 2 (t). For p G [p\(t),<p 2 (t)\, 
s G [t, t] we have v\(s, p) = 0 and v 2 (s, p) < M(s) < £\ by assumption. By Lemma |2.4| , 
for s < t 2 (p) we have dv 2 /dt < —^ and therefore t 2 (p) < min(t + 2 M(t),T). Hence 
(p 2 {t + 2 M(t)) < p\ (f), and the result follows. 
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Lemma 2.9 There is a constant C 2 = C 2 (T, Co) and an increasing function H: [0, T] —> 
M depending on v± and V 2 , satisfying H( 0) = 0 and H(T ) < C 2 , such that if Mft) < £1 
for 0 < t < t, then 


[ |0i(s) - 9 2 (s) I ds < C 2 M( 0) + [ M(s) dH(s) 

•Vo ./o+ 

as Zon# as t + 2Mft ) < r. 

Proof: Using that f Vj^dtp > v- 2 > Gp 2 , from the formula for Off) we obtain that 

101 ft) ~ d 2 ft)\ < C 2 |<^i(t) - <p 2 ft )| + C* / \v\ /Z - U 2 /3 | 

■J 0 

Let ip+(t) = max(<^ 1 (t), (p 2 ff)), then for 95 < <^ + we have 

1 1/3 1/3, ^ \vi - v 2 \ 

-V l< 2 / 3 , 2 / 3 - 

Ui + 1 V 


Note that from Corollary 2J5, it follows that with t+fp) = max(ii(</?), t 2 {}p)) we have 

rMvO 
'o uT 


/•MvO 1 

2 /3 2/3 dt<C(T,C 0 ). 

JO vd 4 - vJ 


By Fubini’s theorem it follows that the function defined by 

!•¥>+(*) 1 


r 

ho(t) = 

■Jo 


2/3 . 2/3 

V + V 




is finite for a.e. t and is integrable with ho ft) dt < C(T, Co). Then we have 

/ — — M(t)hoft) (15) 

ao 

for a.e. t. € [0, T], 


Next, for j = 1 and 2 we invoke Lemma 2.2 with G(f) = —<pjff), fft) = 2 Mft), and 
conclude that as long as t + 2M(t) < r, then 


f <pj(s) - <pj(s + 2M(s)) ds < 2M(0) + f 2M( 
Jo J 0+ 


(s) dHj(s) 


where Hjft) = —<Pjft + 2 Mft)) + pj(2M(0)). Evidently i/j satisfies Hjft) < 1 for all t. 
Putting these estimates together with the result of Lemma 2.8, we find that 


| 6 »i(s) - 0 2 (s)| ds < 4C 2 M(0) + [ M(s) dH(s) 

J 0+ 
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where 


H(t) = 2Cl{H 1 {t) + H 2 (t )) + Cf [ h 0 (s ) ds. 

Jo 


The desired result follows. 


The proof of Proposition 2.6 uses a continuation argument based on the estimates 
above together with the estimate 


M(t) - M{t) < 2 Cl(r - t) 


(16) 


whenever 0 < £ < r < T, which follows from \dv/dt\ < C 3 . Since M is increasing, we 
can find T <T such that T + 2 M{T) = T. With r = t + 2 M(t), inequality ( |l6|) yields 

M(t + 2 M(t)) < M(t)( 1 + 4C 3 ) (17) 


whenever t <T. Now let 


n = {t G [0, T] I M(t + 2 M(t)) < ei}. 


If Af(0) < do := £i/(l + 4C 3 ), then 0 € so 0 is nonempty, and clearly is closed. We 
claim is open in [0,T] if M(0) is sufficiently small. 

~~I to deduce that 


Given any t\ € we can apply Lemmas 2.7 and 


M(t) < Ci(l + C 2 )M(0) + Gi [ M(s ) dtf(s) 

J o+ 


(18) 


for 0 < t < ti. Then Lemma 2.1 implies 


M(t) < C 3 M(0) 


(19) 


for 0 < t < ti, where Cs(T,Co) = exp(C'iC , 2 )C'i(l + C 2 ). Using (p~7| ) we infer that 
M(t\ + 2M(ti)) < C^MiO) with C4 = C 3 ( 1 + 4C 3 ). Provided we assume 


M(0) < <Si := 


l£l_ 
2 (7 4 ’ 


it follows that M(t\ + 2M{t\)) < £ 1 , and since M is continuous, is open in [0, T], 
Consequently we have T G SI, Putting t\ = T, this means we have M(T) < £1 and 
M(T) < C 4 M(0) if M(0) < 5]. This finishes the proof of Proposition |2.6|. 


Proof of Theorem |l.l| : Uniqueness follows immediately from Proposition To prove 
existence for arbitrary vq G rcd{ [0,1]), by Proposition 2.6 and Lemma 2.3 it evidently 
suffices to prove global existence for vq in a dense set of rcd([ 0,1]). Solutions in general 
are constructed by passing to the limit in C([0, T],rcd{[ 0,1])) for every T > 0. 


Lemma 2.10 The set of functions in rcd([0,1]) that take a finite number of values is 
dense in rcd{[ 0,1]). 
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Proof: Let vq £ rcd{[ 0,1]) and let e > 0. Let yj = for j = 0,1, ..and let 

v E (p) = min {yj \ yj > u 0 (</?)} 

for p £ [0,1]. It is easy to see that v £ has a finite number of values, that v £ £ rcd([ 0,1]), 
and ||u e — uo|| < e. This proves the lemma. 


Suppose, then, that vq £ rcd([ 0,1]) takes a finite number of values yo > ■ ■ ■ > Un = 0. 
Then with ipj = \i\i{p \ vq( p) = yj}, we have 0 = po < ■ ■ ■ < Pn < 1 and vo(p) = yj for 
p G [<pj, (fj+i), j = 0,..., N — 1. We start to construct a solution by solving the system 
of ordinary differential equations 


w 


'(t) = Wj(t) 1 / 3 @(t) - 1, j = 0,... ,N- 1, 


( 20 ) 


with 


( 21 ) 


/N—l 

Q(t) = p N / w i( t ) 1/3 ( < Pj+ 1 “ <Pj) 

/ 3=0 

and Wj( 0) = yj, on a maximal interval [0,tjv) in which min Wj(t) > 0. The solution is 
smooth and Wj(t ) > Wj+i(t) by backwards uniqueness for the equation w' = u; 1//3 0 — 1. 
The quantity 


TV—1 

w j( t )( , P3+i - w) 

3=0 

is conserved in time. Without loss of generality we can assume this quantity is 1. 

We can estimate ©(f) < wo(t) 2 ^ 3 so w' 0 < wq and hence wo(t) < e t yo. < oo, then, 
it follows that the smallest component vanishes, i.e., vjn-i (f^) = lin Myt N WN-i(t) = 0. 

For t G [0, tjy) we define v(t, p) = Wj{t) for p G \pj, Pj+i), j = 0,. .., N — 1, and let 
0 = 0. This yields a solution of equations (ED and © for t G [0, t]\r). As t —> tjy from 
below, the limits v(t] y, p) and 6{t] y) exist. The solution can then be re-initialized at time 
tjv with one less component (A^ replaced by N — 1). After some finite number of such 
steps the solution must exist globally. 

Thus, for vq G rcd( [0,1]) with a finite number of values, a global solution exists. 


Theorem 1.1 follows. 


Proof of Theorem L2: The additional restriction imposed in Proposition can be 
removed now by considering convex combinations of initial data. Given T, Co, v\, and V 2 
as stated, let C > 0, <5 > 0 be as given by Proposition | 2 . 6 |, and let Mq = ||ui(0, ■) —(0, Oil- 
Fix an integer n > M.q/5, and for j = 0,1 ,n let 


Xj(p) = (l _ Vl (°’ V) + V2 (°’ 

for p G [0,1]. Then Xj G rcd([ 0,1]), Xj( 0) < Co for all j and ||xj + i — Xj\\ = Mq/ti < 8. 
By the existence theorem there exist corresponding solutions v = Vj to m-m with 
Vj( 0, •) = Xj, and by Proposition ^]6| we have 

sup \\v j+ i(t, ■) - Vj(t, -)|| < C\\x j+ i - Xj\\ = CMo/n. 

0 <t<T 


12 




Since v\ — v 2 = X^=o ('O'+i — Uj), using the triangle inequality we find that 

sup ||ui(f, •) - v 2 {t, -)|| < CM 0 , 

o <t<T 


as desired. 


3 Measure-valued solutions 


Our aim here is to describe a precise correspondence between the solutions v(t,<p ) of 
Theorem IT and measure-valued weak solutions vt of i) , and to show that the metric 
||iq — U 2 II on rcd([0,1]) corresponds to the L°° Wasserstein metric on the space Vo of 
(Borel) probability measures on [0, 00 ) with compact support. Theorem LS then follows 
as a corollary of Theorems |1.1| and |1.2| . 

We begin with a technical lemma on generalized inverses of increasing functions. 


Lemma 3.1 Suppose b > 0 and w: [0,6] —» M is a left continuous increasing function 
with w( 0) =0. Let = w(b ) and define w t; [ 0 , 6 ^] —> R by 


w ] {y) 


supjx | w(x) < y}, 0 < y < b', 

0 , y = 0 . 


Then w ^ is left continuous and increasing, and moreover, 


w'' = w. 


Proof: Clearly is increasing. Given y G (0,6^] and e > 0, put x = w^(y) and 
2 5 = y — w{x — e). Then 5 > 0 and w(x — e) < y — 5, hence x — £ < w\y — 5) < x. It 
follows w t is left continuous. 

To show w ^ = vj, it suffices to show that for 0 < x < b, 

w(x — e) < w'' (x) < w(x + e) 

for all sufficiently small £ > 0. Let y = w''(x) = sup{y | uf(y) < x}. Then for all 
£0 > 0, w^(y + £ 0 ) > x, hence for any small £ > 0 we have w{x — e) < y + e 0 , therefore 
w(x — e) < y. 

For the reverse inequality there are two cases: If x = w'(y) < x then for small £ > 0 
we have 

y < vj(x + e) < w(x + e). ( 22 ) 

Otherwise x > x, and since w is left continuous, x = x. In this case, ( | 22 | ) again holds. 
This finishes the proof. 


If w is continuous and strictly increasing, then uA is the inverse function of w. 

Given a probability measure u with compact support [0,u] C [0, 00 ), we associate the 
distribution function F u : [0, 00 ) —* [0,1] given by 


F„(x) 


z'QO, x)) x > 0 
0 , x= 0 . 


(23) 
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F v is left continuous and increasing, and F u determines v (that is, the values of F v 
determine the values of v on all Borel sets). We associate a decreasing function v = v(v) 
to v via v(x) = Fu( 1 — x) for x £ [0,1]. (Here, Fl is the generalized inverse of the 
restriction of F v to [0,u + 1].) That is, 


( sup{y | F v (y) < 1 - x}, 0 < x < 1, 

\0, x=l. 


(24) 


With the notation Rv{x) = v{\ — x) we have v{y) = R(Fl). The first part of Lemma 

probability measures on [0, oo) with compact support.) 

The inverse map to v is given as follows. If v £ rcd([ 0,1]) we let F = ( Rv )1 on 
[0,w(0)] and put F(x) = 1 for x > u(0). Then F is increasing and left continuous, and 
determines a (Borel) probability measure v for which F = F v . For later use we note that 
for any continuous /: (0, oo) —> R with compact support, we have 


3.1 implies v{y) £ rcd([ 0,1]), thus the map v:Vq —> rcd{[ 0,1]). (Recall Vq is the set of 


/*1 P 1 P OO POO 

/ f(v(x))dx= f(F\x))dx= / f(y) dF{y) = / f(y)dv(y). (25) 
Jo Jo Jo Jo 

This follows from jl], 2.5.18(3)], for example. The identity function y y can be approx¬ 
imated uniformly on compact sets in [0, oo) by such functions /, hence 




U dv{y). 


We let v{v) = v, so i>: rcd{[ 0,1]) —> Vq. Lemma ^]T] implies that we have 


(26) 


Lemma 3.2 v and 0 are inverse maps: v(v(v)) = v for all v £ rcd([ 0,1]), and v(v{v)) = 
v for all v £ Vq. 


We now recall from || that the L p Wasserstein metric can be defined on Vq as follows. 
Given and V2 in Vo, let D(y i, V2) be the set of probability measures y on [0, 00) x [0, 00) 
with marginal distributions z^i and v 2 , that is, for all continuous (: [0, 00) —> R, 


POO POO POO 

/ / c (x)dp(x,y)= C (x)dvi(x) 

Jo Jo Jo 


and 


POO POO POO 

/ / ({y)dy(x,y)= C( y)dv 2 {y )• 

Jo Jo Jo 

If 1 < p < 00 then the L p Wasserstein metric is defined by 

d P {v \,v 2 )= inf / \x-y\ p dy(x,y) 

\/i£D(v 1 ,u 2 ) J 


The L°° Wasserstein metric is defined by 


dooiyi j ^2) = inf |U-ess sup|x — y\. 

p£D(v 1 , 1 / 3 ) 
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The measures /i represent ways to ‘rearrange mass’ from one distribution into the other, 
and the L p Wasserstein metrics measure the least costly way to do this according to the 
notion of cost indicated. 


Lemma 3.3 

have 


Given u\, v 2 E Vo, let v\ = v{y\), v 2 = vfa)- Then for 1 < p < oo we 


r*l \ 1/P 

dp(vi,V2) = ( / \vi(<p) - v 2 (ip)\ p dp ) 


and 


doo(vi,V2) = ||^i - v 2 \\. 

Proof: The assertion for 1 < p < oo follows from 0, see pp 28-30 and Corollary 7.3.6, 
which yields that 

<b{v i ,V2)= Wh (<p) - Fl 2 (ip) \ p d(pj . 

Then Proposition 3 of Q asserts that lim^oo d p {y\, u 2 ) = doo(^i 5 v 2 ). Since v\ and v 2 
are right continuous, it follows 


doolTi, v 2 ) 



\ 1 /p 

v 2 {y)\ p dip) 


ess sup \ vi ( ip ) - v 2 (ip) | = sup \vi (ip) 

[ 0 , 1 ] [ 0 , 1 ] 




Corollary 3.4 Let Vo have the topology induced by d^. Then Vo is complete, and the 
map v'.Vo rcd([ 0,1]) is an isometric isomorphism of complete metric spaces. 

The completeness of Vo with respect to the metric d^ was established in [2] . 

The correspondence between volume orderings v E rcd([0,1]) and volume distribu¬ 
tions v E Vo has been established. Now we seek to show that this correspondence maps 
solutions to weak solutions and vice-versa. 

Proposition 3.5 Let 0 E oo), v E C([0, oo), rcd([0, 1])) be a solution of &■ For 

each t > 0, let v t = 0(v(t,-)). Then u: [0, oo) —> Vo is locally Lipschitz, and v is a weak 
solution of f||). 

Proof: From ([II]) we have that v: [0, oo) — > rcd([0, 1]) is locally Lipschitz, therefore 
v. [0, oo) —► Vo is locally Lipschitz by Corollary |3-4 

Let C : (0, oo) x (0, oo) —► M be smooth with compact support. Then for all <p E [0,1], 
1 1 —> C(t, v(t, ip)) is Lipschitz continuous and we have 

f°° d 

0 = J —({t,v(t,ip))dt 

dtC(t , v(t, ip)) + A (v(t, ip), 0{t))d v C{t , v(t, ip)) dt. 
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Let F(t,-) = (Rv(t,-))^ and let ip(t,-) = 1 — F(t,-). We integrate over ip £ [0,1], use 
Fubini’s theorem, and change variables using (^). We obtain 


0 = - 


(d t ((t,v) + A (v,0(t))d v ((t,v)) d<p(t, v) dt 

roo roo 

/ / (d t C(t,v) + A(v,6(t))d v ((t,v)) dv t {y)dt. 

Jo Jo 


1 0 JO 

r»OQ pOO 


Thus ( 6 , z/) is a weak solution in the sense of Theorem |1.3|. as claimed. 


Proposition 3.6 Suppose that 6 £ L“ c (0, oo) and w. [0, oo) —* Vq is locally Lipschitz, 
and v is a weak solution of &■ Let v(t, •) = v(v t ) for each t > 0. Then v is a solution 

o/ (HP- 

Proof: Given 6 and v as described, the map v: [0, oo) —> rcd([ 0.1]) is locally Lipschitz. 
We consider test functions f of the form 


C (t,v) = C(t)v(v) 


(27) 


where the functions £, p: R —> M are smooth with compact support in (0, oo). Using this 
form together with the fact that (0, u) form a weak solution to (|[), and using the change 
of variables from (25) as previously, we find that 

pOO pi 


0 = 


/ o 


/ p)) + ip))A(v(t, ip), 6(t)) dip dt. 

Jo 


Using Fubini’s theorem and integrating by parts in time, this gives 

n l poo 


0 = 



£(t)f)(v(t,ip))(A(v(t,<p),9(t)) - d t v{t,ip))dtdip. 


(28) 


o Jo 


where fj = rf. This formula is justified since for each ip £ [0,1], v(-,ip) is Lipschitz, 
hence differentiable almost everywhere. We note that since v is bounded on compact 
sets, rj can be chosen to agree on the range of v with an arbitrary smooth function with 
compact support in (0, oo). We do this and drop the tilde. Furthermore, we note that by 
Lebesgue’s dominated convergence theorem, (|2^) remains valid for any £ with compact 
support in (0, oo) that is the bounded pointwise limit £ = lim, woo of a sequence of 
smooth with compact support in (0, oo), and similarly for t). For the moment it will 
suffice to consider £, fj £ (7 C (M + ), the set of continuous functions on (0, oo) with compact 
support. 

For what follows, we take some care regarding joint measurability in (t, ip) and sets 
of measure zero. We fix a representative 0 in the equivalence class 0 , then drop the tilde. 
dtv(t,ip) need not exist at every point, but equation (fis]) also holds if dtv is replaced by 
the upper derivative dtv or the lower derivative d t v, defined by 


dtv(t, ip) 


lim sup 5 h v(t,ip), d t v(t,ip) = lim inf 5 h v(t,ip), 

£ ^°0<|/i|<£ e^0 0<|/i|<e 
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where 


s h u \ v(t + h,p) -v{t,p) 

6 v(t,p) = -^-. 

Since v is locally Lipschitz in t uniformly in p, dtv and d t v are bounded on compact sets, 
and d t v < dtv. 

Lemma 3.7 As maps from (0, oo) x [0,1] —> R, dtv and d t v are Borel measurable. 
Moreover, dtv = d t v almost everywhere in (0, oo) x [0,1]. 

Proof: Since v is continuous in t uniformly in p and is right continuous and decreasing in 
p, v is lower semicontinuous, hence Borel. Suppose 0 < t\ < t 2 < oo, then for 0 < \h\ < t\ 
the map 5 h v is Borel on [ti, t?\ x [0,1]. Let {hj} be a dense sequence in (—1, 0) U (0,1). 
Since the maximum of two Borel functions is Borel and pointwise limits of sequences of 
Borel functions are Borel, and pointwise we have 

sup 6 h v(t,p) = sup S hj v(t,p) = lim max 5 hj v(t, ip), 

0<W<* \hf<e 

by taking e to zero along a sequence it follows that dtv is Borel on [ti,^] x [0,1], hence 
on (0, oo) x [0,1]. A similar argument applies for d t v. 

Now we have that the set Z = {(f, <p) \ (dtv — d ± v)(t, p) > 0} is a Borel set. We know 
that for each p, v(-,p) is differentiable almost everywhere, so ( dtv — d t v)(t,p) = 0 for 
almost every t > 0. Fubini’s theorem now implies that Z has Lebesgue measure zero in 
(0,oo) x [0,1]. 

Returning to (p8|), we can now apply Fubini’s theorem and deduce that for almost 
every t, (dtv — d±v)(t, ! p) = 0 for almost every <p, and with 

J v ( t )= [ v(v{t,<p))(A(v(t,(p),e(t))-d t v(t,(p))dip, 

■Jo 

we have that for any rj € C c (M + ), / 0 °° £(t)J v (t) dt = 0 for all f € C c (M + ). Therefore, 
given i] there is a set C (0, oo) of full measure (meaning the complement has measure 
zero), such that J v {t) = 0 for all f £ 

The set C' c (M + ) is separable, so if we do this for a dense sequence {r] n } we find there 
is a set O C of full measure in (0, oo) such that for t £ O, J Vn (t) = 0 for all n. Since 
any rj £ C' c (M + ) can be approximated uniformly by a subsequence of {ri n }, we infer that: 

Lemma 3.8 There is a set 17 C (0, oo) of full measure, such that for all t £ Q, dtv{t, p) 
exists for almost every p £ [0,1], and J v {t) = 0 for all r) £ C^R" 1 "). 


Lemma 3.9 Let t £ 17, and suppose v{t,x) = v(t,y ) where 0 < x < y < 1. Then 
dtv(t,p ) exists for all p £ (x,y), and is constant on this interval. 

The proof of this lemma is straightforward, using the facts that v(t, •) is decreasing 
for every t, and dtv(t, p) exists for almost every p. 
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Lemma 3.10 Let t £ Ll, and let p(t) = sup{y? | v(t, p) > 0}. Then 


A(v(t, <p),0(t)) - d t v(t, p) = 0 


for almost all ip £ [0, p(t)). 


Proof: We thank B. Kirchheim for the main idea of the following proof. Since t is fixed, 
we suppress indicating dependence on t, and we let g(p) = A (v(t,p),9(t)) — dtv(t,p). 
We know g is measurable and bounded, ip i—>• v{p) is decreasing, so if 0 < y is in the 
range of v, either the pre-image v~ l {y) is a singleton or an interval of nonzero width. 
There can be only a countable set of y of the latter type. Let A be the (countable) set 
of endpoints of such intervals. For x, y £ [0,1] \ A, we know that v(x) = v(y) implies 
d t v{x) = d t v(y), and so g(x) = g(y). 

Given any e > 0, let A s = [0, </?(t))n{a; £ [0,1] | g(x) > e}\A. Then A e is measurable, 
and we claim the measure of A e is zero for any e > 0. Suppose not, so \A e \ = 28 > 0 
for some e > 0. By Lusin’s theorem, there is a compact K C A e such that \K\ > 5 and 
v \k is continuous. Then v(K) is compact and is contained in (0, oo) since v is positive 
at each point of [0, ip(t)). 

Apply Lemma with r)(<p) = g n {p) = max{0,1 — n dist(<^, n(A"))} for n = 1,2,... 
Then r] n has compact support in (0, oo) for n sufficiently large and converges boundedly 
pointwise to the characteristic function Xv(K)- It follows that 


0 = 



Xv(K)(v(ip))g(p) dp. 


Now if v(x) £ v(K), then v(x) = v(y) for some y £ K, and either g(x) = g(y) > e or 
x € A. It follows 

l 

Xv(K){v(p))g(p) dp > e\K\ > 0, 

yielding a contradiction. Hence \A e \ = 0 for any e > 0. A similar argument applies for 
{x | g(x) < — e}, and we then deduce that g{p) = 0 for almost every p £ [0, p(t)). This 
proves the Lemma. 


Note that for t H and p £ (p(t), 1] we have that v(t, p) = 0, and d tv(t, p) = 0. 
Since (t,p) i—>• v(t,p ) is right continuous in p and locally Lipschitz in t uniformly in 
p, the set Q = {(f, p) \ v(t, p) > 0} is open in (0, oo) x [0,1]. Define 

g(t, p) = A (v(t, p),0(t)) - d t v(t, p), 

then g is measurable on (0, oo) x [0,1] and by Lemma [3.10| we have 

roc r 1 

/ / XQ\g\dpdt = 0. 

Jo Jo 

By Fubini’s theorem, we have \Q9 = 0 almost everywhere. Hence there exists a set S 
of full measure in [0,1], such that if p £ S then (xQ9)(t, <p) = 0 for almost every t, and 
t > v(t, p) is differentiable almost everywhere. 
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Lemma 3.11 If p G S and v(t, p) > 0, then v(s, p) > 0 for all s G [0, t\ and 


v(t,p) =v(0,p) + / A (v(s,p),9(s))ds. 


( 29 ) 


Proof: For any t\ G (0 ,t), if v(s,p) > 0 for all s G [ti,t] then since s i—> v(s,p) is 
differentiable and g(s,p ) = 0 almost everywhere in we have 


v(t, p) = v(t u <p) + / A (v(s,p),9(s))ds. 

Jti 

We claim that the set {t\ G [0 ,t) \ v(s,p) > 0 for all s G [ti,t]} has the infimum i* = 0. 
Note that the set is nonempty by the continuity of s i—> v(s, ip). Suppose the infimum t* 
is positive. Then v(t*,tp) = 0 < v(s,ip) for s G We know that 9(s) is bounded for 

s G [0, t], so for some sufficiently small h > 0 it follows that A (v(s,ip),9(s)) < —1/2 for 
t* < s < t* + h. Then we have 

rt* -\-h 

0 < v(t* + h, p) = 0 + / A(v(s,ip),9(s)) ds < —h/2, 

Ju 

a contradiction. Hence our claim holds: v(s,F) > 0 for all s G [0, t], and the formula 
asserted in the Lemma follows. 

Now we can finish the proof of Proposition [O]. Suppose t > 0, tp G [0,1] are arbitrary 
and v(t,ip) > 0. Since v(t, •) is right continuous and decreasing, there exists a sequence 
of numbers p n G S such that p n > p, p n —* p as n —>• oo and v{t, <p n ) > 0. Using Lemma 


3.11 it follows that v(s,<p) > 0 for all s G [0 ,t] and 

v(t, p) - n(0, <p) - [ A (v(s,p) 1 9{s))ds 
■Jo 

= lim (v(t,<p n )-v(0,p n )~ / A(v(s,<p n ),9(s))ds ) = 0 . 

rwoo ^ ,/ 0 J 


This completes the proof of Proposition 3.6. 


The results asserted in Theorem 1.3 now follow directly from Theorems 1.1 and 1.2 


with the help of Propositions |3.5| and 3.6, Corollary |3.4|, and equation (26). 


4 A different conserved quantity 

In the theory of Ostwald ripening, one also encounters an alternative to the condition 
that the total particle volume is conserved in time. If mass in the diffusion field is taken 
into account, one finds that a quantity of the form 


Q = a< 


9(t) + [ v(t,<p)dip 

Jo 


(30) 
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is conserved instead, where a > 0 is a constant. 8(t) need no longer be positive. 

In terms of the theory developed in this paper, the constraint ( |30| ) is simpler to deal 
with that the constraint of constant volume. One has the bound 


6(t) < Q/a, 

and when comparing two solutions of one can use the arguments of Lemma |2.7| and 
replace the use of Lemma |2.9| by the simpler estimate 


\0i(t) - d 2 (t)\ < a 1 \\vi(t, •) -v 2 (t, -)|| +a 1 |Qi-Q 2 |- 


-ii 


(31) 


From the standard Gronwall’s inequality, one easily deduces the a priori estimate asserted 
in the following result. The existence and uniqueness proofs go the same as in section 3. 


Theorem 4.1 Let vq G rcd([0, 1]), Q G R. Then there exists a unique function v G 
C([0, oo), rcd([0,1 ])) such that, with 6{t) determined by §3(\ ), we have 

v{t, ip) = v 0 (ip) + f (v(s, p) 1/3 0(s) - 1 ) ds 
Jo 

whenever v(t, p) > 0. 

Given T > 0, Co > 0, there exists a positive constant C such that, given two solutions 
as above which also satisfy Q 2 ) < Co, then 

sup |Mt, ■) - v 2 (t, -)|| < C (||ui(0, •) - u 2 (0, -)|| + \Qi - Q 2 \). 

0 <t<T 


Using the correspondence v(t, •) v t = 0(v(t, ■)) and its inverse as in section 3, from 
Propositions 3.5 and |3.6| one may deduce directly the following corollary of Theorem I 


Theorem 4.2 Given vq G Vq, Q G R, there exists a unique map t 1— > vt that is locally 
Lipschitz from [0, 00) into Vo such that, with d(t) determined by the relation 



vdv t (v), 


we have 

d t C(t,v) + A(v,0(t))d v C(t,v) du t {v) dt = 0 

for all smooth f: (0, 00) x (0, 00) —► R with compact support. 

Given any T > 0, Co > 0, there exists C > 0 such that, if two such weak solutions 
, z/ 2 ) are given, which satisfy max(Qi,(5 2 ) < Co, then 



sup doo(^ t (1) , 4 2) ) < C 
0 <t<T v 


^ 2) ) + \Ql-Q2\). 
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